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Abstract. This paper examines the mode identification 
and presents pulsation models for FG Vir, for which 24 fre- 
quencies have been detected. Histograms of the frequency 
spacings show peaks which are identified with adjacent 
radial orders and rotational splitting. 

Pulsational £ values are deduced for eight modes by 
comparing the observed photometric phase lags between 
v and y variations with calculated values. The dominant 
pulsation mode at 12.72 c/d can be identified with 1=1, 
while the 12.15 c/d mode is the radial fundamental. These 
results are in agreement with identifications published by 
Viskum et al. (1998). 

Based on the observational mode identifications and 
the Hipparcos distance, new models were computed with 
the constraint that the mode at 12.15 c/d is the radial fun- 
damental mode. It is shown that with standard opacities, 
models in the appropriate T ff, logL and log g ranges can- 
not reproduce the identification in the literature of 23.40 
c/d as the third radial overtone. However, we show that 
observationally an t = 1 (rather than radial) identification 
is equally probable. 

A large number of pulsation models were computed 
for FG Vir. A comparison between the observed frequen- 
cies and mode identifications and pulsation models leads 
to a mean density of p/p~Q — 0.156 ± 0.002 depending on 
the opacity and chemical composition choice and on the 
possible overshooting from the convective core. The mod- 
els also correctly predict the observed region of instability 
between 9 and 34 c/d. 

The effect of rotational coupling on the pulsation fre- 
quencies is estimated. 

Key words: 5 Set - Stars: oscillations - Stars: evolution 
- Stars: interiors - Stars: individual: FG Vir - Stars: in- 
dividual: HD 106384 
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1. Introduction 

FG Vir (=HD 106384) is a 5 Scuti variable near the end 
of its main-sequence evolution. 435 hours of photometric 
measurements by the Delta Scuti Network determined 24 
statistically significant frequencies from 9.20 to 34.12 c/d 
(106 to 395 /iHz). Details of this campaign as well as ref- 
erences to earlier measurements and results can be found 
in Breger et al. (1998). The large number of detected pul- 
sation modes makes this star an excellent candidate for 
asteroseismological investigations. This requires the iden- 
tifications of the observed pulsation frequencies with spe- 
cific pulsation modes. While the problem is rather com- 
plex, considerable progress has been achieved, as shown 
by Breger et al. (1995), Guzik et al. (1998) and Viskum 
et al. (1998). 

The pulsation mode indentification from observed fre- 
quencies requires accurate determinations of the basic pa- 
rameters of the star. From the available uvbyfi photome- 
try, Mantegazza et al. (1994) derived T cS = 7500 K and 
log g = 3.95. A correction for a misprint in the literature 
leads to a correction of log g to 3.9. We can now improve 
these values further by including the accurate Hipparcos 
parallax which predicts My = 1.95 ± 0.13 mag. This is 
slightly fainter than the value of 1.71 ±0.25 mag predicted 
by uvbyP photometry. This leads to a corresponding shift 
in log g to 4.00. These values are in exact agreement with 
those derived by Viskum et al. (1998). We estimate the 
uncertainties to be ±100 K in temperature and ±0.1 in 
log 9- 

The values of the pulsation constants Q can be esti- 
mated from the following empiric equation: 

logQ, = -6.456 + logPi + 0.51og 3 + 0.1 M ho} + logT cff . 

The constant, -6.456, in the above formula is based on 
solar values of Afboi = 4.75 mag, B.C. — —0.08 mag, T e s = 
5770 K and log g = 4.44. If the Q values are calculated 
from uvbyf3 photometry, the observational uncertainties 
in observing these parameters lead to an uncertainty in Q 
of about 18%. 
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Table 1. Pulsation frequencies of FG Vir 



Frequency 1995 V amplitude Q value 

c/d fiRz mmag days 



Statistically significant frequencies 



/l 


12.716 


147.2 


21.1 


.0323 


h 


24.228 


280.4 


4.5 


.0170 


h 


23.403 


270.9 


4.1 


.0176 


h 


21.052 


243.7 


3.7 


.0195 


h 


19.868 


230.0 


3.5 


.0207 


h 


12.154 


140.7 


3.5 


.0338 


h 


9.656 


111.8 


3.4 


.0426 


h 


9.199 


106.5 


3.1 


.0447 


h 


19.228 


222.5 


1.5 


.0214 


ho 


20.288 


234.8 


1.3 


.0203 


hi 


24.200 


280.1 


1.3 


.0170 


fl2 


16.074 


186.0 


1.0 


.0256 


fl3 


34.119 


394.9 


1.0 


.0121 


/l4 


21.232 


245.7 


1.0 


.0194 


/IB 


11.110 


128.6 


0.9 


.0370 


fw = 2/i 


25.432 


294.4 


0.9 


.0162 


hr 


33.056 


382.6 


0.6 


.0124 


fia 


21.551 


249.4 


0.8 


.0191 


fl9 


28.140 


325.7 


0.6 


.0146 


ho 


11.195 


129.6 


0.7 


.0367 


hi 


24.354 


281.9 


0.6 


.0169 


/22 


11.870 


137.4 


0.4 


.0346 


/23 = fl + /7 


22.372 


258.9 


0.5 


.0184 


/24 = /s - /l 


10.687 


123.7 


0.5 


.0385 


Probable frequencies 








/25 


25.37 


293.7 


0.4 


.0162 


/26 


25.18 


291.4 


0.4 


.0163 


/27 


29.50 


341.4 


0.4 


.0139 


/28 


18.16 


210.2 


0.4 


.0226 


/29 


19.65 


227.4 


0.4 


.0209 


/30 


31.92 


369.4 


0.4 


.0129 


/si 


20.83 


241.1 


0.4 


.0197 


/32 


12.79 


148.1 


0.4 


.0322 



The corresponding Q values are shown in Table 1. 

2. Regularities of frequency spacing 

The values of the observed frequencies and regularities in 
their patterns can be an excellent intital tool for mode 
identifications, if enough frequencies arc excited and de- 
tected. For high-order, low-degree p-mode pulsation, the 
different radial orders show uniform frequency spacing, 
with a mode of order n and of degree I being shifted 
from the corresponding mode (n, £+ 1) by half of the fre- 
quency difference between the (n,£) and (n + 1,£) modes 
(Vandakurov 1967). In S Scuti stars, the excited pulsa- 
tion modes are of low order (n up to 7), so that the 
asymptotic relations do not apply exactly. Nevertheless, 
they also show some regularities. Additionally g-modes 
invade the p-mode region and decrease the spacing in a 



small frequency region of about two radial orders. This 
effect, known as avoided crossing (Osaki 1975, Aizcnman, 
Smcyers & Weigert 1977), complicates the theoretical fre- 
quency spectra, but can provide information about the 
stellar interior (Dzicmbowski & Pamyatnykh 1991). More- 
over, stellar rotation splits multiplets and this splitting is 
non-symmetric, if second-order effects of rotation and ef- 
fects of rotational mode coupling are taken into account 
(Dziembowski & Goode 1992, Soufi et al. 1998, Pamyat- 
nykh et al. 1998). Nevertheless, the spacing of adjacent 
radial orders as well as the rotational splitting is still reg- 
ular enough to be detectable, if complete multiplets are 
excited and identified. We will demonstrate this by using 
a pulsation model of a 1.85 M Q star with T e s = 7515 K, 
log g = 3.99, and V ro t — 45km/s. This model will be re- 
ferred to as model 1. The parameters for the model were 
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Fig. 1. Histograms of frequency spacing between all specified pulsation modes. Top left: The diagram demonstrates 
that for high orders the patterns of frequency spacing clearly show adjacent radial orders (~ 4 c/d) and the effects of 
rotational splitting, which is extremely asymmetric even at V ro t — 45 km/s. Top right: The frequency spacing predicted 
from model 1 for i = 1 in the observed frequency range of 11 - 35 c/d. Note that the patterns from adjacent orders 
and rotational splitting arc still present. Bottom panels: Observed frequency spacings in the observed range from 11 
to 35 c/d. Although these are a mixture of t = 0, 1 and 2 modes, the effects of adjacent radial orders and a small peak 
in the range of rotational splitting can be seen. To demonstrate that the results are not sensitive to which observed 
frequencies are included, two different choices (see text) are shown 



not chosen at random, but can be regarded as an estimate 
for FG Vir. 

To investigate the period regularities, Winget et al. 
(1991) have successfully applied the method of the Fourier 
transform of the period spacing to the star PG 1159+035. 
This method requires coherence over a large frequency 
range. Handler ct al. (1997) also found frequency regular- 
ities from Fourier transformations of the frequency spec- 
trum of the unevolved <5 Scuti star XX Pyx. Since strict 
equidistant frequency or period spacing is not expected for 
FG Vir, the method is not optimal for this 6 Scuti star. 
Instead, we use a method which does not require such a 
coherence: an examination of a histogram of the observed 
frequency differences between all detected frequencies. In 
such a diagram, regularities in the frequency spacing of 
adjacent radial orders of modes with the same degree, £, 



should show up as a peak. Furthermore, modes of different 
degree are shifted in frequency relatively to each other, but 
would still have similar patterns and, therefore, contribute 
to the peaks in the histogram. 

The frequency spacing is examined in Fig. 1 with 
both the theoretically predicted and observed spacings. 
Pulsation models show a typical frequency spacing of 
A/ « 4 c/d for adjacent radial orders of p-modes, inde- 
pendent of the degree of the modes. The leftmost peaks 
in the top panels of Fig. 1 are caused by rotationally split 
multiplets. A similar diagram for £ = 2 (not plotted sep- 
arately) does not show such strong peaks in the expected 
region. The reason is that both the presence of g-modes 
in addition to the p-modes and non-equidistant rotational 
splitting significantly disturb the regularity in the distri- 
bution of quadrupole mode frequencies (see Fig. 7 below.) 
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As a result, the combined pattern of frequency spacings 
for all I = — 2 modes becomes much less clear. More- 
over, due to the fact that only low-order oscillations are 
present in this frequency range, there is no additional peak 
at Af « 2 c/d as might be expected from the asymptotic 
spacing between p-modes of adjacent degrees (see Fig. 7 
for more details). 

Next, we turn to the observed frequency spacing for the 
24 certain and 8 probable frequency detections of FG Vir 
(Tabic 1). The most cautious approach would be to use 
the 24 certain frequencies with a few exceptions: the 2/i 
term at 25.4 c/d (reflecting the departure from a pure si- 
nusoidal light curve shape of /i), the two combination fre- 
quencies (the pulsation models cannot yet predict which 
combinations and resonances are excited), and the two 
low-frequency modes for which the p-mode character can 
definitely be excluded from the assumption that fe is the 
radial fundamental mode. To show that the agreement be- 
tween the theoretically predicted and observed frequency 
spacing is not based on the choice of frequencies, the anal- 
ysis was repeated by including the 8 additional 'probable' 
modes listed in Table 1. 

To conclude, the theoretical and observed frequency 
spacings agree quite well. In particular, for frequency dif- 
ferences in the 0-5 c/d range, two features near 3.9 and 
0.8 c/d stand out, suggesting an identification with the 
spacing of successive radial orders and rotational splitting, 
respectively. 

3. Pulsation mode identifications from 
photometric phase differences 

The relative phase difference between the temperature and 
radius variations of a pulsating star leads to an observ- 
able phase difference between the light curves at differ- 
ent wavelengths. The sizes of these phases differences de- 
pend not only on the properties of the star, but also on 
the type of pulsation mode. The observed phase differ- 
ence can then be used for mode typing. This was already 
pointed out by Watson (1988). Garrido et al. (1990) pre- 
sented detailed calculations and predictions for S Scuti 
stars. They find that measurements through different fil- 
ters of the Stromgren uvby system provide discrimation 
between radial and low-order nonradial pulsation, i.e. help 
determine the I value 

We have chosen the v and y filters to provide a rela- 
tively large baseline in wavelength. The u filter was not 
used by us because of the very large potential for system- 
atic observational errors. These details of the measure- 
ments can be found in Breger et al. (1998). The phase 



T = 7500 K, Logg = 4.00 



1 It is necessary to note that rotational mode coupling may 
enlarge the overlapping between modes of different £ values in 
the amplitude-phase diagrams, as it was discussed by Pamyat- 
nykh et al. (1998): for example, a quadrupole mode coupled by 
rotation with the closest radial mode may be shifted in such a 
diagram towards the region occupied by dipole modes. 
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Fig. 2. Diagnostic diagram to determine I values of FG 
Vir from Stromgren v and y colors. The axes represent 
amplitude ratios and phase differences. Measurements are 
shown by crosses with error bars, while the four-sided 
loops represent the models (see text). The three panels 
represent the pulsation modes with different values of the 
pulsation constant, Q 



differences were determined in the following manner: The 
values of the 24 known and well-determined frequencies 
were optimized by making a common solution of the avail- 
able y data from 1992 - 1996, while allowing for the am- 
plitude variability of /a. As discussed in Breger et al., all 
CCD measurements were given a weight of 0.19. With 
these optimized frequencies, for the year 1995 the best 
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Frequency Phase differences in degrees Puisation degree, I 



c/d 


fiRz 


<f>v 


- 4>y 


Spectroscopy 


Photometry 






1995 


1995/6 


Viskum et al. (1998) 


Present 


fi 12.716 


147.2 


-1.0 ± 0.2 


-1.3 ± 0.2 


1 


1 


/ 2 24.228 


280.4 


-3.0 ± 1.1 


-3.8 ± 1.1 


1 


1, 2 


fs 23.403 


270.9 


-0.6 ± 1.2 


-1.1 ± 1.2 





0, 1 


/ 4 21.052 


243.7 


-5.7 ± 1.4 


-7.1 ± 1.5 


2 


2 


fs 19.868 


230.0 


-4.9 ± 1.4 


-5.5 ± 1.5 


2 


2 


/ 6 12.154 


140.7 


+6.8 ± 1.4 


+5.5 ± 1.4 








f 7 9.656 


111.8 


-2.2 ± 1.4 


-2.7 ± 1.4 


2 


1, 2 


fs 9.199 


106.5 


-4.4 ± 1.6 


-7.6 ± 1.7 


2 


2 


Number of hours y/v 


412/292 


494/374 








0.5 - 

, l , l , l , l , l , l , l , l , 

-10 -8 -6 -4 -2 2 4 6 8 

Phase difference, 4> v -<]> y , in degrees 

Fig. 3. Comparison of the equivalent-width and photo- 
metric methods to determine I values. Radial pulsation 
(£=0) can be found in the lower right, 1=1 in the middle, 
while £=2 is found near the top left. The diagram shows 
that the two methods are in agreement, but also demon- 
strates that the some of the agreement may be accidental 
once the error bars are taken into consideration 

amplitudes and phases were calculated from the available 
412 hours of y and 292 hours of v data. Separate trial 
solutions indicate that the resulting phase differences are 
relatively insensitive to the weights adopted. For the year 
1996, an additional 82 hours of uvby photometry are avail- 
able (Viskum et al. 1998). The data were combined with 
the larger data set from 1995 while allowing for variable 
amplitudes of / 3 .We note that the calculated uncertain- 
ties of the phase differences are not reduced by including 
the additional data: the reason is that the 1995 data have 
smaller deviations, e. g. 4 vs. 6 mmag in v. 

The resulting phase differences are shown in Table 2. 
The phase errors (in degrees) were estimated from the 
formula a((f)) = 180/-7T • a(m)/a, where a is the amplitude 
and a(m) is the uncertainty of each data point (average 
deviation per point from the fit). 



We can now compare the observed phase differences 
with theoretical modelling in order to determine the I val- 
ues. The ATLAS9 models of Kurucz (1993) were used to 
construct a model atmosphere for FG Vir. Garrido et al. 
(1990) presented calculations using values of $ T ranging 
from 90° to 135°. For FG Vir, Viskum (1997) determined 
a smaller range, viz. <J> T = 126° ± 20°. This allowed us to 
rchnc the calculations, although the results are very sim- 
ilar. Another required constant, the deviation from adia- 
baticity, R, has been changed slightly from the value used 
by Garrido et al. (1990). Values of 0.20 (instead of 0.25) 
to 1.00 were used. This change was indicated by mea- 
surements of high-amplitude S Scuti stars. The theoretical 
predictions are shown in Fig. 2 together with the observa- 
tions. The importance of considering the dependence on 
the pulsation constant, Q, can be seen for Q = 0.02, where 
one can even find negative values of the phase difference 
for radial modes, although the separation between radial 
and nonradial modes is always maintained. 

Our best mode identifications based on Stromgrcn 
photometry are shown in Table 2. We obtained five unam- 
biguous I values, while for three further modes we cannot 
distinguish between two adjacent I values. The frequency 
/ 6 , shown in the middle panel, is situated to the right of 
the I = by 1.9 a. We note that the deviation is caused 
by only one subset of data (the CCD measurements from 
Siding Spring Observatory, see Stankov et al. 1998), with- 
out which a phase shift of +3.7 ° is found. Irrespective of 
which of the two values for f§ is accepted, an identification 
with radial pulsation is consistent within the statistical 
uncertainties. 

We can now compare the results from the photometric 
method with those derived from a promising new tech- 
nique of examining the equivalent width variations of se- 
lected lines. Bedding et al. (1996) have shown that for low 
degree pulsation, the ^-values of pulsation modes can be 
inferred from simultaneous observations of several selected 
absorption lines combined with simultaneous photomet- 
ric observations. Viskum et al. (1998) have applied this 
method to the star FG Vir. In particular, the equivalent- 
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width changes of the Ha and H/3 lines turned out to be 
good discriminators. In their paper, i identifications have 
been presented for the eight dominant modes. 

We note that on the observational side the photomet- 
ric and spectroscopic methods are independent. However, 
both methods rely on similar model-atmosphere calcula- 
tions, so that they cannot be considered to be completely 
independent of each other. 

The agreement between the photometric and spectro- 
scopic mode determinations is remarkable. It appears pru- 
dent to examine the comparison of the results of the two 
methods in more detail, especially with consideration of 
the (unavoidable) observational uncertainties. In order to 
compare independent parameters with each other, we pick 
the amplitude ratio of A(Ha)/A(FeI) given by Viskum et 
al. (1998). The comparison is shown in Fig. 3, where the 
numbers next to the points refer to the frequency num- 
bering in Table 1. The figure shows that some of ex- 
cellent agreement may be accidental once the observa- 
tional uncertainties are considered. Nevertheless, the vi- 
ability of both methods to determine i values has been 
demonstrated and for at least six modes the I values 
have been observationally determined. These determina- 
tions now need to be used as input for pulsation models. 

4. Pulsation models for FG Vir 

Since the initial discovery of multiperiodicity of FG Vir, 
several studies attempted to fit the observed and theoreti- 
cal frequency spectra of the star, viz. Breger et al. (1995), 
Guzik, Tcmplcton & Bradley (1998), and Viskum et al. 
(1998). We will now calculate new models utilizing the 
newly discovered pulsation frequencies and mode identifi- 
cations. 

Method of computation 

To compute models of FG Vir we used a standard stel- 
lar evolution code which was developed in its main parts 
by B. Paczyhski, M. Kozlowski and R. Sienkiewicz (pri- 
vate communication). The same code was used in our re- 
cent studies of period changes in <5 Scuti stars (Breger & 
Pamyatnykh 1998) and in a seismological study of XX Pyx 
(Pamyatnykh et al. 1998). These two papers include de- 
tailed descriptions of the model computations, so that the 
present decription can be brief. For the opacities, we used 
the latest version of the OPAL or the OP tables (Iglesias & 
Rogers 1996 and Seaton 1996, respectively) supplemented 
with the low-temperature data of Alexander & Ferguson 
(1994). In all computations the OPAL equation of state 
was used (Rogers et al. 1996). 

The computations were performed starting with chem- 
ically uniform models on the ZAMS, assuming typical 
Population I values of hydrogen abundance, X, and heavy 
element abundance, Z. The initial heavy element mixture 
of Grevesse & Noels (1993) was adopted. 



In some models, a possibility of overshooting from the 
convective core was taken into account. The overshooting 
distance, rf vor, was chosen to be 0.2 H p , where H v is the 
local pressure scale height at the edge of the convective 
core. Examples of evolutionary tracks for <5 Scuti mod- 
els computed with and without overshooting are given in 
Breger & Pamyatnykh (1998). 

In the stellar envelope the standard mixing-length the- 
ory of convection with the mixing-length parameter a = 
1.0 or 2.0 was used. As we will see below, the choice of 
the mixing-length parameter a has only a small effect on 
our models, because they are too hot to have an effective 
energy transfer by convection in the stellar envelope. 

In all computations we assumed uniform (solid-body) 
stellar rotation and conservation of global angular momen- 
tum during evolution from the ZAMS. These assumptions 
were chosen due to their simplicity. The influence of ro- 
tation on the evolutionary tracks of <5 Scuti models was 
demonstrated by Breger & Pamyatnykh (1998). We stud- 
ied models of FG Vir with equatorial rotational velocities 
from, approximately, 30 to 90 km/s (on the ZAMS, the 
values are 5-10 km/s higher). This range is consistent with 
the values of v sin i = 21 ± 1 km/s and i = 31° ± 5° found 
by Mantegazza et al. (1994) and an equatorial velocity of 
33 ± 2 km/s obtained by Viskum et al. (1998). At such 
low rotational velocities, the evolutionary tracks are lo- 
cated very close to those for non-rotating stellar models. 
The main effect of rotation to be considered is the split- 
ting of multiplets in the oscillation frequency spectra. This 
splitting is non-symmetric even for slowly rotating stars, 
if second-order effects are included. 

The linear nonadiabatic analysis of low-degree oscilla- 
tions (£ < 4) was performed using the code developed by 
Dzicmbowski (1977). In the modern version of the code, 
effects of slow stellar rotation on oscillation frequencies 
are taken into account up to second order in the rota- 
tional velocity (Dziembowski & Goode 1992, Soufi et al. 
1998). 

4-2. Model constraints using oscillation data 

The models for FG Vir were constructed with the ob- 
served mode Jq (12.154 c/d) being identified with the ra- 
dial fundamental mode (— F) (see Section 3). Note that 
this determines the mean density of all possible models 
of FG Vir: with the pulsation constant of about 0.032 - 
0.034 days, which is typical for <5 Scuti variables, we obtain 
p/f>Q « 0.15-0.17. A considerably more accurate value of 
the density will be obtained later in this section. 

We started with the construction of evolutionary 
tracks of 1.75 - 1.95 M Q models for initial abundances 
X = 0.70 and Z = 0.02 and using OPAL opacities. No 
overshooting from the convective core was allowed. The 
initial equatorial rotational velocity on the ZAMS was 
chosen to be 50 km/s. With our assumption of conser- 
vation of global angular momentum, the equatorial rota- 
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tional velocity is decreasing during the MS-evolution from 
50 km/s at the ZAMS to about 40-41 km/s at the TAMS 
(Terminal-Age-Main-Sequence). The evolutionary tracks 
are shown in Fig. 4 together with the range in effective 
temperature and gravity of FG Vir (see Introduction) de- 
rived from photometric calibrations (T c g = 7500 ± 100 K 
and log g = 4.00 ± 0.10). This range requires MS models 
and constrains the mass of the models to 1.75-1.95 M Q . 
The position of the models, whose F-mode has a value 
of 12.154 c/d is also shown. This further constrains the 
mass to 1.815-1.875 M . We stress that this strong seis- 
mological mass constraint depends on an accurate effective 
temperature determination. 

The identification of / 6 = F agrees well with the grav- 
ity estimate for FG Vir. This provides no additional con- 
straints on the stellar mass since the lines of constant fre- 
quencies are approximately parallel to those of constant 
gravity, as shown in the lower panel of Fig. 4. 

For a given family of stellar models, the radial fun- 
damental pulsation constant, Q, is constant with a quite 
high accuracy due to the homologous structure of mod- 
els of different masses. For M = 1.80 - 1.90 M & in the 
range log T cS = 3.869 - 3.881 (T cff = 7400 - 7600 K), the 
pulsation constant Q = 0.0326 with relative accuracy of 
about 0.2 %. The accuracy will be still higher by approxi- 
mately one order of magnitude if we consider only models 
based on / 6 = F. For these models we determine a mean 
density of p/p© = 0.1570 ± 0.0001. However, such an ex- 
tremely high accuracy is based on a fixed choice of input 
physics: stellar opacity, initial chemical composition, rota- 
tional velocity and parameters of convection. We will see 
in the next subsection that changing these parameters re- 
sults in significantly larger spread of mean density for FG 
Vir models. 

The strict constraint on mass is demonstrated in Fig. 5, 
where the changes of radial and dipole frequencies during 
MS evolution are shown. This agreement is an indepen- 
dent qualitative argument in favour of the proposed mod- 
els for FG Vir. Another important result is the good agree- 
ment of the predicted frequency range of unstable modes 
with the observed frequency range of 9-34 c/d. An addi- 
tional test shows that should we identify fe with the first 
radial overtone instead of the F-mode, we cannot achieve 
agreement between the theoretical and observed frequency 
ranges: in the corresponding models of w 2.0M Q the in- 
stability occurs in the frequency range of 8-30 c/d. The 
tendency in models of higher mass to shift the instability 
range to lower frequencies can be also seen in Fig. 5. There 
is an even stronger argument against these higher-mass 
models: their luminosities are too high (log(L/L ) <~ 1.3) 
to be consistent with both the photometric calibrations 
and the Hipparcos parallax. 

A number of gravity modes must be excited in FG Vir, 
if the assumption of f§ = F is true, because the two low- 
est frequencies are more than 25% lower than the F-mode. 
Moreover, during the MS-evolution the frequencies of low- 




3.9 3.85 
log T*, 

Fig. 4. Evolutionary tracks of 1.75-1.95 Mq standard 
models. The equatorial rotational velocity on the ZAMS 
was chosen to be 50 km/s. On the TAMS (at the turn- 
off points to the left), equatorial velocities are of about 
40-41 km/s. Dashed lines show effective temperature and 
log g ranges of FG Vir from photometric calibrations. The 
thin solid line connects models whose radial fundamental 
mode frequency is 12.154 c/d. The rotational velocity of 
the models along this line is about 45 km/s 



order g-modes increase and approach consecutively the 
frequencies of low-order p-modes resulting in mode inter- 
actions and avoided crossings (see Unno et al. 1989 and 
references therein) . The frequency spectrum is much more 
complicated than in the case of pure p- or g-modes, as 
shown in Fig. 5 for dipole modes. The avoided crossing 
phenomenon takes place approximately in the middle of 
the observed frequency interval. Therefore, most of the ex- 
cited modes at these and at lower frequencies are of mixed 
character: they behave like p-modes in the envelope and 
like g-modes in the interior. In the 1.85 M Q model with 
log T c ff = 3.876, modes g\, p 2 and p 3 are of mixed char- 
acter. The frequencies of modes at avoided crossing are 
sensitive to the structure of the deep stellar interior. Con- 
sequently, the detection of these modes is important for 
testing convective overshooting theories (Dziewbowski & 
Pamyatnykh 1991). 
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Fig. 5. Main-sequence evolution of low-order frequency spectra of radial and dipole oscillations of stellar models with 
masses 1.80, 1.85 and 1.90 M Q . In each panel, lefmost and rightmost points correspond to the ZAMS and to the TAMS 
models, respectively. Large filled circles denote unstable modes. For simplicity, for I = 1 only axisymmetric (m = 0) 
components of the dipole multiplets are shown. Rectangular boxes mark the observational frequency and effective 
temperature range of FG Vir. The vertical line in each panel denotes a model whose radial fundamental frequency 
(F-mode) fits the observational frequency — 12.154 c/d. Only models with masses 1.815-1.875 M Q fit the allowed 
temperature range 



Avoided crossings for quadrupole modes in the mod- 
els of FG Vir occur close to the upper border of the ob- 
served frequency interval and also close to the F-mode. 
This means that most of I = 2 p-modes in the interval 
already interacted with gravity modes and are of mixed 
character. 



4-3. Effects of different input parameters on the FG Vir 
models 

The 1.85 M Q model for FG Vir, which was discussed in 
the previous subsection, will be referred to as the standard 
or reference model with the input parameters: X = 0.70, 
Z = 0.02, d OVCI = 0, a = 1.0, Kot.zAMS = 50 km/s and 
OPAL opacities. To examine the effects of varying input 
parameters on the predicted frequency spectrum, all these 
and the stellar mass were varied, under the condition that 
F = ft. 
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Table 3. Parameters of FG Vir models with F = f$. The symbols have their usual meaning (see text). For the opacity, 
k, the OPAL, OP or arficially modified OPAL data were used, pi/p^ is the ratio of frequency of radial fundamental 
mode, /(pi), to that of third overtone, f(jpA) 



Model 


M/Mq 


X 


Z 


logT cfi 


logL 


R/Rq 


log g 


Vrot 


a 


Clover 




p/Pe 


Pl/P4 


1 


1.85 


0.70 


0.02 


3.8760 


1.1690 


2.274 


3.988 


45 


1.0 


0.0 


OPAL 


0.1571 


0.5236 


2 


1.82 


0.70 


0.02 


3.8701 


1.1406 


2.261 


3.985 


45 


1.0 


0.0 


OPAL 


0.1571 


0.5233 


3 


1.85 


0.70 


0.02 


3.8761 


1.1696 


2.274 


3.989 


31 


1.0 


0.0 


OPAL 


0.1560 


0.5231 


4 


1.85 


0.70 


0.02 


3.8756 


1.1676 


2.274 


3.983 


67 


1.0 


0.0 


OPAL 


0.1570 


0.5248 


5 


1.85 


0.70 


0.02 


3.8753 


1.1656 


2.272 


3.976 


90 


1.0 


0.0 


OPAL 


0.1563 


0.5266 


6 


1.85 


0.70 


0.02 


3.8760 


1.1691 


2.274 


3.988 


45 


2.0 


0.0 


OPAL 


0.1570 


0.5236 


7 


1.85 


0.70 


0.02 


3.8796 


1.1837 


2.275 


3.987 


45 


1.0 


0.2 


OPAL 


0.1558 


0.5231 


8 


1.85 


0.65 


0.03 


3.8734 


1.1562 


2.267 


3.990 


45 


1.0 


0.0 


OPAL 


0.1584 


0.5227 


9 


2.00 


0.70 


0.03 


3.8748 


1.1844 


2.327 


4.001 


46 


1.0 


0.0 


OPAL 


0.1574 


0.5231 


10 


1.72 


0.70 


0.02 


3.8754 


1.1507 


2.233 


3.972 


45 


1.0 


0.0 


OP 


0.1542 


0.5262 


11 


1.95 


0.70 


0.02 


3.8712 


1.1600 


2.301 


4.000 


45 


1.0 


0.0 


mod.OPAL 


0.1588 


0.5204 


12 


1.95 


0.70 


0.02 


3.8746 


1.1738 


2.301 


4.002 


32 


1.0 


0.2 


mod.OPAL 


0.1597 


0.5195 



The changes introduced by using different opacities or 
non-standard chemical composition were mainly compen- 
sated by changes in mass, in order to fulfill the only iden- 
tification we use. 

The main characteristics of twelve models of that se- 
ries are given in Table 3. Model 2 differs from model 1 
(our reference model) in mass; models 3, 4, 5 - in rota- 
tional velocity; model 6 versus 1 will demonstrate effect 
of changing the mixing-length parameter a; model 7 ver- 
sus 1 will show effect of the overshooting; models 8 and 9 
have non-standard chemical composition; finally, models 
10, 11 and 12 differ from model 1 in opacity (additionally, 
overshooting is taken into account in the model 12). 

Note the significantly larger spread in stellar mass be- 
tween different models (1.72-2.00 Mq) than for the mass 
interval of 1.815-1.875 Mq in the case of the standard 
choice of input parameters as was discussed in the previous 
subsection. The same is true for the mean density range: 
in Table 3 it varies between p/p® = 0.1542 and 0.1597 
(or between 0.1558 and 0.1584 when using only OPAL 
opacities). This spread is at least one order of magnitude 
larger than for the standard input data. Nevertheless, this 
seismic estimate of the mean density, which is based both 
on the well determined effective temperature and the one 
identification we are using, provides a strong constraint 
on possible FG Vir models [| 

We note that besides quite different stellar masses of 
1.7 - 2.0 Mq (see Table 3) the evolutionary tracks for 
all 12 models in their MS-part lie well inside the region of 
1.80 - 1.90 Mq of the standard set. Including a luminosity 
estimation log(L/L ) w 1.1 — 1.2 from trigonometric par- 
allax determined by Hipparcos, all MS model tracks pass 

2 For the multiperiodic 8 Scuti-type star XX Pyxidis, for ex- 
ample, there is no observational information about mode iden- 
tification. Therefore it was necessary to consider a large num- 
ber of models with very different mean densities (Pamyatnykh 
et al. 1998). 



the error box, as well as the error box in the log <?Tog T e g- 
diagram. On the contrary, none of the post-MS models fits 
such a combination of parameters. 

4-4- The problem of the radial frequency ratio 

Viskum et al. (1998) identified f^ as a radial mode. We 
note here that the phase-difference method presented ear- 
lier in this paper allows both £ — and 1 for /3, i. e. radial 
as well as nonradial pulsation. We will now examine the 
radial hypothesis. In the last column of Table 3 the ratio 
of frequencies of the radial fundamental mode, p\ , and of 
the third overtone, p4, is given {f{pi)/f(pi))- For mod- 
els 1-10 these values are close to, but not equal to the 
observed ratio of fe/ fs = 0.5193, independent of which 
parameter was changed. This can also be seen in Fig. 6, 
where the ratio f{pi)/f(Pi), is plotted against f{p\) for a 
wide range of parameters of S Scuti star models. There are 
well-defined monotonic variations of this ratio with chang- 
ing mass, effective temperature or chemical composition, 
but the observed ratio disagrees with all these results. 

The only exception is the model 12 with artificial opac- 
ities, which was constructed in the following way. For FG 
Vir models, the frequency ratio is most sensitive to the 
choice of opacities as it can be seen from comparison of 
models 10 and 1 (OP versus OPAL opacities). Physically, 
OP data differ from those of OPAL by underestimation 
of collective effects in stellar plasma, therefore OP opac- 
ities are systematically lower than OPAL in deep stellar 
interiors. For FG Vir models, this difference in opacity 
is about 20% at temperatures above 10 6 K. In the en- 
velope, at lower temperatures, OP opacity varies slightly 
more monotonously along radius than does OPAL opac- 
ity: some dips are slightly shallower and some bumps are 
more flat. The differences do not exceed 8%: for example, 
at a temperature of 14 000 K the OP opacity is 4% smaller 
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8 10 12 14 16 

radial fundamental frequency, f(pi), [c/d] 



Fig. 6. Frequency ratio of the radial fundamental mode to the third overtone for a wide range of parameters of 5 Scuti 
star models and of some FG Vir models from Table 3 (asterisks) . The large filled circle corresponds to the observed 
frequency ratio, fe/f3, of 0.5193. Only the models with artificially modified opacities (such as model 12 of Table 3) 
can fit the observed ratio 



and at a temperature of 300 000 K it is 7.5% higher than 
the OPAL opacity. 

Using the fact that the difference in frequency ratios 
between model 1 (OPAL) and model 10 (OP) is compa- 
rable with the difference between model 1 and the ob- 
servations (but these differences are of opposite sign, see 
Fig. 7), we performed a very simple numerical experiment: 
we artificially scaled OPAL opacities with a factor, which 
is the ratio of OPAL to OP data. More clearly, we used 
^modified = kopal • [kopal/kop]- Models 11 and 12 were 
constructed using K mo dified- For model 12 we additionally 
set d vor to 0.2 H p and lowered the rotational velocity. 
This model fits the observed frequency ratio very nicely as 
demonstrated in Fig. 7. However, this agreement should 
not be construed as an indicator for a new revision of 
atomic physics data on opacity, since the mode identifi- 
cation from Viskum et al. (1998) may not be unique due 
to the size of the error bars. Moreover, we cannot exclude 
additional effects like nonlinear mode interaction or rota- 
tional mode coupling, which may influence the frequency 
spectrum. In the last section the problem of rotational 
mode coupling will be briefly discussed. The observed vari- 
ability of the amplitude of mode fs is another argument 
in favour of possible nonlinear mode interaction. 

Viskum et al. (1998) were able to interpret the ob- 
served frequency ratio fe/fs = 0.5193 as the radial fre- 
quency ratio f (pi) / f (pi) ■ They did not construct full evo- 
lutionary models but scaled a model of 2.2M Q , which was 
selected to match the observed frequency ratio with the 
radial frequency ratio f(p\)/f{pi)- Using the homology 
argument, they estimated the mean density of the true 
FG Vir model by multiplying the mean density of the 



2.2Mq model by the square of the ratio / O bs//modoi- In 
such a way an agreement between observed frequencies 
/6, /3 and a pair of radial modes of the scaled model was 
achieved by definition. The estimated gravity, luminos- 
ity and distance of the scaled model were found to be in 
good agreement with the photometric and the spectro- 
scopic data and with the Hipparcos parallax. The authors 
noted that the high precision of their asteroseismic den- 
sity estimate (p/p~Q — 0.1645 ± 0.0005) is based on a fixed 
(solar) metallicity for FG Vir. Indeed, with our standard 
choice of chemical abundances and opacity and assuming 
the fit F = fe we estimated the mean density of the FG Vir 
model with even five times higher accuracy (see section 
4.2), but possible variations of the global parameters re- 
sult in at least one order of magnitude worse precision of 
this estimate. 

4-5. Theoretical frequency spectra versus observations 

Frequency spectra of radial, dipole and quadrupole modes 
for all 12 models from Table 3 are shown in Fig. 7. The 
effects of different choices of input parameters can be es- 
timated by comparison of the results for different models. 
We discuss here both general properties and some pecu- 
liarities of these frequency spectra. 

For nonradial oscillations, evolutionary overlapping of 
frequency intervals of g- and p-modes (see Fig. 5) re- 
sults in avoided crossings, which disturb the approxi- 
mately equidistant frequency spacing between acoustic 
multiplcts. Gravity and mixed modes are very sensitive 
to the interior structure as can be seen for models of dif- 
ferent chemical composition, different opacities and for 
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Fig. 7. Frequency spectra of radial, dipole and quadrupole oscillations of various FG Vir models. Axisymmetric modes 
(m = 0) arc marked by enlarged circles. Model numbers (see Table 3 for parameters) together with some model 
indicators are given to the right of the panels. Vertical solid and dashed lines correspond to observed frequencies - 
statistically significant and probable, respectively. Numbers above some observed frequencies give identifications for 
the degree of the modes (£) based on multicolor photometry data and on the results by Viskum et al. (1998) 



models with and without overshooting. On the contrary, 
the change of a (model 6 versus model 1) has a negligi- 
ble influence on the frequency spectrum due to ineffective 
convection in the relatively hot envelope of FG Vir. 

Rotation splits nonradial multiplets and strongly com- 
plicates the frequency spectra. Except for models with 
slowest rotation, we observe a forest of quadrupole modes 
in the low-frequency part of the interval, with overlapping 
components of the different multiplets. The common prop- 



erty of the spectra is a large asymmetry of the rotational 
splitting, which is caused by the second-order effects of 
rotation (Dziembowski & Goode 1992). The asymmetry 
is higher the higher the order of the p-modes is. 

It is not trivial to select a model reproducing the ob- 
served frequencies exactly. Simple attempts to minimize 
frequency differences (O-C) by a combined variation of 
input parameters of the reference model fail due to strong 
and non-linear sensitivity of gravity modes to interior 
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structure. At the same time, this strong sensitivity may 
help to fit some chosen frequencies without changing the 
rest of the spectrum, (cf. models 1 and 7, for example). 
It is obvious from Fig. 7 that generally it is much easier 
to fit a low-frequency mode than a high-frequency mode, 
because the spectrum is more dense at lower frequencies. 
For some of the observed modes there is no satisfactory 
£ = — 2 solution: see, for example, the group f^, /14, /ig 
around 21 c/d or fn at 33 c/d. Besides geometric cancel- 
lation, there are no objections to identify frequencies in 
the gaps with modes of degree £=3 and £=A. Note that 
even for t = — 2 the number of unstable modes is a 
few times larger than the observed one: in the observa- 
tional frequency interval there are 6-7 radial modes, 24 
dipole modes and 50-55 quadrupole modes. Therefore, a 
presently unknown mode selection mechanism must exist. 

Note that most of the models presented in Fig. 7 show 
a good fit of the dominant observed mode /1 (12.716 c/d) 
with £ = 1 mode of m = —1 or 0. This is in agreement 
with the mode identification from photometric phase dif- 
ferences. On the contrary, it is quite difficult to achieve a 
similar fit with a dipole mode for the observed mode f% 
(23.403 c/d). In Table 4 we present some results to quan- 
tify the fitting of 21 modes (corresponding to the observed 
frequencies /1 through / 22 , omitting f 1& ) for models 1 (ref- 
erence model), 3 (low rotation), 8 (X = 0.65, Z = 0.03) 
and 12 (artificially modified opacity + overshooting). In 
the cases of close observed frequencies (for example, /1, 
hi, /21 around 24 c/d, or / 4 , f 14 , f 18 around 21 c/d) 
we give a few possible identifications for each frequency. 
Moreover, if there is no £ = — 2 mode for a given ob- 
served frequency, we show the closest mode of £ = 3 or 
4: the frequency spectrum of these modes is so dense that 
practically everywhere in the observed interval a fitting 
within 0.1 c/d is possible. 

Model 12 seems to be the best-fitting model in our 
series. Note the excellent agreements in both frequency 
and ^-identifications for all ten dominant frequencies. The 
mean difference (O-C) is about 0.04 c/d for these modes. 
The fit for most of the remaining frequencies is also good, 
once £ = 3 £ = 4 modes are considered. Possible dis- 
crepances at low frequencies do not appear serious to us 
due to the strong sensitivity of g-modes to model param- 
eters. Because of the effects of avoided crossing some of 
the quadrupole modes of higher order are also quite sen- 
sitive to small variations of parameters (cf. £ = 2 modes 
for models 1 and 2 at approximately 32 c/d). 

In Fig. 8 we compare the range of observed frequencies 
and of excited modes explicitly for the reference model 
as well as three other models. FG Vir is located in the 
HR-diagram in the middle of the instability strip (see, 
for example, Breger & Pamyatnykh 1998). Consequently, 
the driving of oscillations is effective in a wide frequency 
region which extends over 7 radial overtones. The inde- 
pendence of the driving efficiency on the mode degree, £, 
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Fig. 8. Normalized growth-rates, rj, of low-degree oscilla- 
tion modes of some FG Vir models from Table 4. Only ax- 
isymmetric modes (m = 0) are shown. Positive values cor- 
respond to unstable modes. Vertical lines mark observed 
frequencies from Table 1, probable frequencies /2s to / 32 
are shorter. The longest line at 12.716 c/d corresponds to 
the mode /1 with the highest amplitude 

is typical for the oscillations excited by the opacity mech- 
anism. 

We also note that the observed frequency spectrum is 
divided clearly into two groups: at 9-13 c/d and at 19- 
25 c/d. However, the present models do not predict any 
instability gap in frequency, as can be seen in Fig. 8. 

4-6. Rotational mode coupling 

An additional factor which influences the frequencies of 
a rotating star is the coupling between modes with close 
frequencies whose azimuthal orders, to, are the same and 
whose degrees, £, are also the same or differ by 2. The 
effect was described and discussed in detail by Dziem- 
bowski & Goode (1992) and by Soufi et al. (1998). The 
frequency distance between two near-degenerate modes 
increases when coupling is taken into account. The sig- 
nificance of this rotational frequency perturbation was 
demonstrated by Pamyatnykh et al. (1998) in applica- 
tion to XX Pyx. It was shown that at rotational velocity 
of about 90 km/s the frequency shifts of coupled radial 
and quadrupole (or dipole and octupole) overtones achieve 
0.1-0.2 c/d. Therefore, in particular, a significant change 
of radial frequency ratios may be expected. 
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Table 4. Best fits for some FG Vir models 



Observations 






Model 1 








Model 3 








Model 8 








Model 12 




N 


/ [c/d] 


id. 


£ 


m 


/ [c/d] 


|0-C| 


£ 


m 


/ [c/d] 


|0-C| 


£ 


m 


/ [c/d] 


|0-C| 


£ 


m 


/ [c/d] 


|0-C| 


1 


12.716 


1 


1 


-1 


12.821 


.105 


1 


-1 


12.714 


.002 


1 


-1 


12.814 


.098 


1 


-1 


12.719 


.003 








2 


-1 


12.770 


.054 


























2 


24.228 


2/1 


1 





24.436 


.208 


1 


1 


24.083 


.145 


1 





24.453 


.225 


1 


1 


24.250 


.022 








3 


-2 


24.331 


.103 


3 


-2 


24.128 


.100 


3 


-2 


24.381 


.153 


3 


-2 


24.308 


.080 


3 


23.403 











23.226 


.177 








23.237 


.166 








23.247 


.156 








23.402 


.001 








3 


1 


23.285 


.118 










3 


1 


23.317 


.086 










4 


21.052 


2 


1 


-1 


20.828 


.224 


1 


-1 


20.758 


.294 


1 


-1 


20.667 


.385 


1 


-1 


20.739 


.313 








2 


2 


21.767 


.715 


2 


2 


22.005 


.953 


2 


2 


21.723 


.671 


2 


2 


22.106 


.054 








3 


-2 


21.079 


.027 


3 


-3 


21.092 


.040 


3 


-2 


20.961 


.091 


3 


-3 


21.024 


.028 








4 


-3 


21.048 


.004 


4 


1 


21.091 


.039 


4 


-3 


21.037 


.015 


4 


-4 


21.085 


.033 


5 


19.868 


2 


2 


-2 


19.996 


.128 


2 


-2 


19.800 


.068 


1 


1 


19.893 


.025 


2 


-2 


19.774 


.094 
























2 


-2 


19.899 


.031 










6 


12.154 











12.161 


.007 








12.156 


.002 








12.152 


.002 








12.157 


.003 


7 


9.656 


1/2 


1 


-1 


9.609 


.047 


2 


-2 


9.705 


.049 


1 


-1 


9.250 


.406 


2 


2 


9.656 


.000 








2 


-1 


9.556 


.100 










2 


1 


9.815 


.159 
































3 


-3 


9.661 


.005 










8 


9.199 


2 


1 


1 


9.228 


.029 


2 





9.267 


.068 


1 


-1 


9.250 


.051 


1 


-1 


9.143 


.056 








2 





9.242 


.043 










2 


-2 


8.940 


.259 










9 


19.228 


(0) 


1 





19.204 


.024 


1 





19.217 


.011 


2 





19.265 


.037 


2 





19.326 


.098 


10 


20.288 


(1) 


1 


1 


20.162 


.126 


1 


1 


20.293 


.005 


1 





20.385 


.097 


1 


1 


20.206 


.082 








3 





20.333 


.045 


























11 


24.200 


- 


1 





24.436 


.236 


1 


1 


24.083 


.117 


1 





24.453 


.253 


1 


1 


24.250 


.050 








3 


-2 


24.331 


.131 


3 


-2 


24.128 


.072 


4 


1 


24.074 


.126 


3 


-2 


24.308 


.108 


12 


16.074 


- 


1 





16.191 


.117 


1 





16.170 


.096 


1 





16.124 


.050 


1 





16.211 


.137 








2 


2 


15.946 


.128 


3 


-1 


16.079 


.005 


3 


-1 


16.016 


.058 


2 


2 


16.123 


.049 








4 





16.060 


.014 


4 





16.082 


.008 


4 


-1 


16.010 


.064 


4 


2 


16.114 


.040 


13 


34.119 




2 





34.231 


.112 


2 





34.221 


.102 


2 


1 


33.916 


.203 


2 


2 


33.970 


.149 
























3 


-2 


34.172 


.053 


3 





34.019 


.100 


14 


21.232 




1 


-1 


20.828 


.404 


1 


-1 


20.758 


.474 


1 


-1 


20.667 


.565 


1 


-1 


20.739 


.493 








3 


-2 


21.079 


.153 


3 


-3 


21.092 


.140 


2 


2 


21.723 


.491 


3 


-3 


21.024 


.208 








4 


-4 


21.368 


.136 


4 





21.367 


.135 


4 





21.167 


.065 


4 





21.147 


.085 


15 


11.110 




2 


-1 


10.975 


.135 


2 


-2 


11.099 


.011 


2 


2 


11.423 


..313 


2 


2 


11.649 


.539 








3 


2 


11.099 


.011 


3 


-2 


11.074 


.036 


3 





10.942 


.168 


3 


-1 


11.024 


.086 








4 





11.092 


.018 


4 





11.085 


.025 


4 


4 


10.983 


.127 


4 


2 


11.059 


.051 


17 


33.056 




2 


2 


33.316 


.260 


1 


-1 


32.609 


.447 


1 


-1 


32.785 


.271 


1 


-1 


32.959 


.097 








4 


-1 


33.020 


.036 


4 


-2 


33.152 


.096 


3 


1 


33.092 


.036 


4 





33.032 


.024 


18 


21.551 




2 


2 


21.767 


.216 


2 


2 


22.005 


.454 


2 


2 


21.723 


.172 


2 


2 


22.106 


.555 








3 


-3 


21.433 


.118 


4 


-1 


21.642 


.091 


4 


-1 


21.566 


.015 


4 


-1 


21.421 


.130 


19 


28.140 




1 


1 


27.959 


.181 


1 


1 


28.125 


.015 


2 


-1 


28.247 


.107 


2 


1 


28.062 


.078 








4 


1 


28.205 


.065 


3 


-2 


28.133 


.007 


3 


-1 


28.101 


.039 


3 


-1 


28.163 


.023 


20 


11.195 




2 


-2 


11.277 


.082 


2 


-2 


11.099 


.096 


2 


2 


11.423 


.228 


2 


2 


11.649 


.454 








3 


-2 


11.281 


.086 


3 


-3 


11.309 


.114 


3 


-1 


11.289 


.094 


3 


3 


11.223 


.028 








4 


2 


11.231 


.036 


4 


3 


11.181 


.014 


4 


-2 


11.285 


.090 


4 


4 


11.175 


.020 


21 


24.354 




1 





24.436 


.082 


1 





24.413 


.059 


1 





24.453 


.099 


1 


1 


24.250 


.104 








3 


-2 


24.331 


.023 


3 


-3 


24.349 


.005 


3 


-2 


24.381 


.027 


3 


-2 


24.308 


.046 








4 





24.398 


.044 










4 





24.457 


.103 


4 


1 


24.339 


.015 


22 


11.870 




2 


1 


12.029 


.159 


2 


2 


11.859 


.011 


2 


1 


11.818 


.052 


2 


1 


11.915 


.045 








4 


4 


11.846 


.024 










3 


-3 


11.969 


.099 











We estimated this effect in some of our FG Vir models mode, f{pi), and of the third overtone, f(p4), are changed 

and found it to be unimportant at rotational velocities of due to coupling with closest axisymmctric quadrupolc 

about and less than 45 km/s. For example, for the ref- modes by -0.0035 c/d and 0.0091 c/d, respectively. The 

erence model 1, the frequencies of the radial fundamental effect is much stronger for more rapidly rotating models: 
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for model 5 with V TO t — 90 km/s, the radial fundamen- 
tal mode and third overtone are shifted by -0.0436 c/d 
and 0.1545 c/d, respectively, which results in the change 
of the frequency ratio from 0.5266 to 0.5212. As another 
example, we were able to reproduce the observed ratio 
f&l h = 0.5193 as the radial frequency ratio for a model 
with initial abundances X = 0.65, Z = 0.03 and with V m t 
= 91-92 km/s. However, the rapid rotation seems to be 
rather improbable for FG Vir with v sinz = 21 ± lkm/s, 
as it was discussed by Viskum et al. (1998). 

Note that the coupling effect is higher for higher over- 
tones: for model 5, the shift of the radial sixth overtone 
frequency (34.306 c/d) is 0.520 c/d due to interaction with 
the closest axisymmetric quadrupole mode of 34.219 c/d 
which is shifted by the same quantity 0.520 c/d in the 
opposite direction. Moreover, t = 2 modes are affected 
by coupling with £ = 4 modes, and so on. We conclude 
that for rapidly rotating models it is necessary to take 
the rotational coupling into account in attempts to fit the 
observed frequency spectrum with the theoretical one. 

The rotational coupling results also in a mutual con- 
tamination of amplitudes of spherical harmonic compo- 
nents of interacting modes (Soufi et al 1998, see examples 
in Table 4 of Pamyatnykh ct al. 1998). This adversely 
affects the mode discrimination by means of multicolor 
photometry (see the footnote in section 3) and should 
influence spectroscopic determinations as well. However, 
this effect was found to be important only in models with 
more rapid rotation than found for FG Vir. A more de- 
tailed discussion of the rotational mode coupling prob- 
lem in connection with the interpretation of the observed 
multifrcquency spectrum will be given by Dziembowski & 
Goupil (1998). 
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